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ABSTRACT 

We develop a non-linear relativistic model of the Shapley Concentration (SC) and its 
environs, including the Great Attractor (GA) and the Local Group (LG). We take the 
Shapley concentration as a major attractive centre, and we use the Lemaitre-Tolman 
model. We constrain our model with present day observations, plus the requirement 
that it have a physically reasonable evolution from small perturbations at last scat- 
tering. We investigate possible mass and velocity distributions, and we find that the 
peculiar velocity maximum near the SC is ~ 800 km/s inwards, the density between 
GA and SC must be about ~ 0.9 times background, the mass of the GA is probably 
4 — 6 x fO 15 Mq, the SC's contribution to the LG motion is negligible, and the value 
of the cosmological constant is not significant on this scale. 
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1 INTRODUCTION 

The Great Attractor and the Shapley Concentration are the 
largest and most important cluster concentrations in the lo- 
cal Universe. The Great Attractor is mostly hidden behind 
our Milky Way (Woudt & Kraan-Korteweg 2001). The first 
measurements of peculiar velocities of the galaxies in this 
region showed that the galaxies move towards the Great 
Attractor (Lynden-Bell et al. 1988). Although the motion 
towards the Great Attractor is significant, there is no evi- 
dence for any backside infall onto the Great Attractor. This 
suggests that the galaxy flow in this region is just a part of a 
larger flow, caused by some more massive attractive center. 
Behind the Great Attractor there is the Shapley Concen- 
tration. The Shapley Concentration is not in the Zone of 
Avoidance and hence is more suitable for observations. 

Fig- LH presents the schematic position of the Local 
Group (LG), the Great Attractor (GA) and the Shapley 
Concentration (SC) in a coordinate system with the X di- 
rection in the galactic plane pointing towards the GA, and 
the Z axis pointing to galactic north. Figs. 2-4 depict 
the Shapley concentration with surrounding galaxies, in per- 
pendicular slices 1600 km/s thick and 32000 km/s square. 
Fig. shows the galaxy distribution in a constant X slice 
through the SC, Fig. [3] in a slice of constant Y, and Fig. 21 
a slice of constant Z. The galaxy data were taken from the 
NASA/IPAC Extragalactic Database . As one can see, there 
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are still regions in which there are hardly any observations, 
while many strips covered by very detailed galaxy surveys 
are evident. Furthermore, the largest amount of data are 
from our immediate neighborhood, deceptively suggesting 
that we are living in highly overdense region. Fig. [3] is the 
only slice that includes the local group, and it shows clearly 
how well measured the nearby regions are, how patchy the 
coverage is further out, and how the north is better sur- 
veyed than the south. Unfortunately, individual surveys do 
not cover all the region of interest to us. Thus there is sig- 
nificant uncerainty in the galaxy distribution around the 
SC, but we feel there is enough data to make an interesting 
model. However Figs. 2-4 do show that the SC is a high 
density region, surround by regions of lower density, with 
a number of overdensities in various directions further out, 
such as the GA and the Hercules supercluster (visible at the 
top of Fig. 2). 

This paper provides a new approach to modeling these 
structures. The innovation is that we model them as a one 
connected system, and we intend the models to be gener- 
ally consistent with the astronomical data about the Shap- 
ley Concentration and the surrounding regions. We create 
a relativistic non-linear model of the Great Attractor and 
the Shapley Concentration using the Lemaitre-Tolman (LT) 
metric, which assumes spherical symmetry about some cen- 
tre. We view the SC as a major attractive centre, so we place 
the spherical origin at the SC, and we consider a sphere that 
includes the LG. 

Although Fig. |3 shows the core of the SC is not of 
spherical shape, this is at least partly an artefact of the 
non-uniform observational coverage in this region. On large 
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distance scales however, as can be seen in Figs. [5] and [I] this 
assumption is not a bad first approximation to the galaxy 
data. We further justify the assumption of an overdense 
sphere around the SC at a radius of about 14000 km/s, 
with a less dense region between, as follows. We take the 
geocentric NED data, calculate the number of galaxies in a 
sequence of redshift/velocity bins, and fit it with a smooth 
polynomial. The deviation of this polynomial from a con- 
stant value represents the decreasing fraction of galaxies 
that are detected with distance. We also find the average 
galaxy density in the north galactic cap (10° < b < 90°), in 
the zone of avoidance ( — 10° < b < 10°), and in the south 
galactic cap (—90° < b < —10°). Next we calculate the num- 
bers of galaxies in a set of redshift/velocity shells centered 
on the SC, but we normalise the galaxy numbers at each 
point using the polynomial and the average densities in the 
3 galactic lattitude zones. The resulting histogram is shown 
in Fig. |S] The density maximum at the radius of the GA is 
quite distinct. Also, in practice the galaxies are in a num- 
ber of concentrations distributed around the sphere, such at 
the GA itself, so the actual density in each concentration 
is higher. We intend that our model will be most accurate 
for regions around the SC, the GA and the LG, and it will 
naturally be considerably less accurate for regions further 
out than the SC, because of the scarcity of data there. 

We aim to address the following questions: What is the 
mass and density distribution around the Shapley Concen- 
tration and the Great Attractor? What is the peculiar veloc- 
ity and density distribution in the region between the GA 
and SC? Is there a backside infall towards the GA? Does 
the SC have a significant impact on the motion of the LG? 
Finally we are interested in investigating the influence of the 
cosmological constant on the dynamics of the local Universe. 
Additional constraints on our model are derived from con- 
sidering its evolution. Since it is based on an exact solution 
of the Einstein field equations, we are able to trace the evo- 
lution backward (and forwards) in time (without any limit 
on the size of the density fluctuations) , and we require the 
model to evolve from small initial fluctuations, preferably 
without forming shell crossings. 

The structure of this paper is as follows: in Sec. [5] 
we present the observational data about this region, in 
Sec. [3] the Lemaitre-Tolman model is presented, and Sec. 
0| presents the results of our calculations. 



2 OBSERVATIONAL CONSTRAINTS 
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Figure 1. The schematic distribution of the Local Group (LG), 
the Great Attractor (GA) and the Shapley Concentration (SC) 
in space. The coordinates X, Y, Z are as described in the text. 
The three surfaces correspond to Figs. 2-4 
















3 

x 10* 



Figure 2. The galaxy distribution in a 1600 km/s thick slice 
of constant X passing through the SC. The SC is at the "O", 
galactic north is upwards, and our galaxy is about 10400 km/s 
behind the central "+". See Fig. for the slice position. 



Direct astronomical observations of galaxies provide us with 
their positions on the sky and their redshifts. With vari- 
ous methods of distance estimation (such as Tully-Fisher 
for spiral galaxies, D — a for elliptical galaxies or the SBF 
method) it is possible to convert redshifts into the peculiar 
velocities of the measured galaxies. Once this is done, the 
masses of clusters can be estimated. Masses can also be es- 
timated from the flux of HI emission or X-rays. Below we 
summarise the astronomical data about the Shapley Con- 
centration and the surrounding regions, with which we wish 
our model to be consistent. 

(i) Peculiar velocities: 

• The motion of the Local Group with respect to the 



CMB is 627 km/s (Tonry et al. 2000). The contributions 
to this motion from the gravitational influence of the GA 
and the SC have been estimated by many astronomers. 
Hoffman et al (2001) estimated that motion of the LG 
towards the GA is in the range 100 to 200 km/s. Smith et 
al. (2000) estimated that 50 ± 10% of the LG's motion is 
generated by the GA and the SC. The X-ray observations 
(Kocevski, Mullis & Ebeling 2004) imply that the SC has 
a significant impact on the LG's motion. 

• Tonry et al. (2000) conclude that the SC does not 
have a significant influence on the LG's motion, but their 
results show that the LG's motion towards the GA is in 
the range 200 to 300 km s~\ and the infall velocity at a 
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Figure 3. The galaxy distribution in a 1600 km/s thick slice of 
constant Y passing through the SC. The SC is at the "O", our 
galaxy is at the central "+" , and up is galactic north. See Fig. 
for the slice position. 
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Figure 4. The galaxy distribution in a 1600 km/s thick slice of 
constant Z passing through the SC. The SC is at the "O", our 
galaxy is about 7200 km/s behind the central "+", and the X 
direction corresponds to galactic longitude £ = 320° . See Fig. 
for the slice position. 
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Figure 5. The number density of galaxies in shells centred on 
the Shapley concentration (arbitrary units), normalised to com- 
pensate for the general decrease in the galaxy detection rate with 
distance from earth, as well as the differences in survey coverage 
in the north and south galactic regions and the zone of avoidance. 
The SC is on the left, and the LG at the right. 



distance of 17 Mpc from the center of the GA is ~ 900 
km s _1 . 

• Observations of clusters in the GA and the SC region 
imply a smaller variation in the peculiar velocity profile 
(Lucey, Radburn-Smith & Hudson 2005). According to 
their research, most clusters associated with the GA, have 
positive peculiar velocities (away from us) of ~ 400 km 
s _1 ; clusters that are immediately beyond the GA have 
small peculiar velocities and do not show any evidence of 
backside infall; peculiar velocities of clusters in the region 
between GA and SC are of very small positive amplitude. 

• Observations of a few Snla beyond 100 Mpc show 
high positive values of the peculiar velocity, around the 
1500 km s _1 , however the error bars are large as well 
(Lucey, Radburn-Smith & Hudson 2005). 

The data above imply two possible profiles for the peculiar 
velocity along the line through the SC — a high amplitude 
velocity profile (HVP) and a low amplitue velocity profile 
(LVP). Fig. HJ presents both these profiles. 

(ii) Mass distribution: 

• We expect the current mass distribution to have 
evolved from small initial fluctuations that existed at last 
scattering. Therefore, given the dust equation of state, the 
total mass of our system today is the same as the mass 
of the same comoving sphere at last scattering, in other 
words the average density over the whole region is the 
same as the density of the background cosmology. 

• Bardelli et al. (2000) estimated that the SC has mass 
5.25 x 10 15 M Q within a sphere of radius 14.03 Mpc, and 
mass 8.25 x 10 15 M Q within 19.6 Mpc, based on the ob- 
served excess of galaxies. 

• X-ray observations (Filippis, Schindler & Erben 
2005) suggest lower mass values; specifically 2.35 x 
10 15 M within a radius of 13.6 Mpc, and 2.88 x 10 15 M 
within 16.7. However, these estimates are based on flux 
observations, and so the contribution from the matter be- 
tween the clusters and in the outer regions of the clusters 
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Figure 6. The main present day peculiar velocity profiles used 
in our models. Distance is measured from the LG, with the SC 
on the right, for easy comparison with published estimates. 
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Figure 7. The main present day density distributions for the 
Shapley Concentration and surrounding regions, as used in our 
models. Here distance is measured from the SC, with the radius 
of the LG on the right. 



is neglected. Therefore these estimates can be treated only 
as lower limits. 

• The mass of the GA is hard to determine, since it 
is in the Zone of Avoidance, and estimates vary from 
2 x 10 15 M (Staveley-Smith 2000) to 3.75 x 10 16 M 
(Lynden-Bell et al. 1988). 

• The density contrast in the GA lies between 8 ~ 0.73 
(Tonry et al. 2000) and o w 1.1 (Kolatt, Dekel & Lahav 
1995). 

• Although the region between the GA and the SC is 
not well known, it can be seen from Figs. [^J— |^|that the SC 
is surrounded by low density regions. The density contrast 
in one of them — the Bootes void — is in the range 5 w 
-0.84 to -0.66 (Day et al. 1990). 

Astronomical data suggest several density profiles. Let us 
consider two of them. Namely, HM with high mass values 
for the SC and the GA, and LM with lower masses for both. 
These profiles are presented in Fig. [7] Both these profiles 
include a low density region between the SC and the GA. 



3 THE LEMAITRE-TOLMAN MODEL 

The Lemaitre-Tolman (LT) model is the spherically sym- 
metric solution of Einstein's equations with a dust source 
(Lemaitre 1933, Tolman 1934). In comoving and syn- 
chronous coordinates, the metric is: 



R' 2 (r,t) 
l + 2E(r) 



dr -ir(i,r)dfi 



(1) 



where dfi 2 = d0 2 +sin 2 6&<j) 2 , R' = dR/dr, and E(r) is an ar- 
bitrary function of r. Because of the signature (+, — , — , — ), 
this function must obey E(r) > — |. 

The Einstein equations can be reduced to the following: 



upc 



2M' 
R 2 R' 



\A 2 = 2E(r) 



^ + ±Atf, 



(2) 
(3) 



where R = dR/dt, M(r) is another arbitrary function, k = 
8nG/c 4 , and we take V(t,r) — R(t,r) as the velocity of a 
worldline. 

Equation J3J can be solved by a simple integration: 



dR 



2E(r) + m. + I AT? 2 



= c{t-t B (r)), 



(4) 



where the arbitrary function £s(r) appears as a function of 
integration. This means that the Big Bang is not simulta- 
neous for all worldlines, as in the Friedmann models, but 
occurs at different times for different coordinate distances 
from the origin. 

Thus, the evolution of the LT model is determined 
by three arbitrary functions: E(r), M(r) and tsir). The 
metric and all the formulae are covariant under arbitrary 
(but monotonic) coordinate transformations of the form 
r = f(r'), and by use of such a transformation any one of 
these functions can be given many different forms. Therefore 
the physical initial data for the evolution of the LT model 
consist of two arbitrary functions. 

Apart from the Big Bang singularity that occurs when 
R = along any given constant r worldline (i.e. at t = 
ts{r)), there are also shell crossings at R' — 0, M 1 7^ where 
the density becomes infinite. Shell crossings can be avoided 
by setting the initial conditions appropriately (Hellaby & 
Lake 1985). 

3.1 The model setup 

We model a spherical inhomogeneous region centred on the 
SC, with a comoving LT radius r max that includes the local 
group. We require that this region matches smoothly, in the 
Darmois sense, to a reasonable "background" Robertson- 
Walker (RW) model at this outer radius. 

The chosen background model is the RW model with 
density 



p b = 0.27 x p c 



0.27 



8ttG' 



(5) 
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where the Hubble constant is -ffo = 72 km _ Mpc - . 
When models with a cosmological constant are considered, 
the value of A corresponds to £2a = 0.73, where Qa = 
(l/3)(c 2 A/tf 2 ). 

For the matching of the LT interior to the RW back- 
ground, if So and po are the present day scale factor and 
density, r;, is the standard RW radial coordinate, and rb m 
its value at the matching surface, then 2E(r rnax ) = — fer 2 m , 
M(r max ) = Kp (S r bm ) 3 /6, and t B {r m ax) = (i.e. we set 
t = at the bang in the homogeneous background.) Thus, 
the relative density inside the LT region is p(t, r)/pb(t), and 
the peculiar velocity of a worldline relative to the SC is 
V(t,r) — Vb{t,r) = R(t,r) — rSo(t). However, for easy com- 
parison with published peculiar velocity profiles, we plot 
the peculiar velocity relative to the earth, which, between 
the LG and the SC, is the negative of the foregoing, i.e. 
r§ (t) - R(t,r). 

For most of the following models, the sphere radius cor- 
responds to the distance between the Shapley Concentration 
and the Local Group, r max = J'lg (about 14500 km/s ~ 
205 Mpc today), but in our last investigation it is slightly 
larger. 

Methods for extracting the arbitrary functions that 
specify a LT model from given data such as density or ve- 
locity profiles have been given in Krasihski & Hellaby 2002, 
2004a, 2004b; Bolejko, Krasihski & Hellaby 2005; Hellaby 
& Krasihski 2006. Specifically, a given initial density profile 
p(ti,r) at time t\ plus a given final density profile p(t2,r) at 
time ti fully determine a particular LT model. In that work 
the coordinate choice r — M was made. Similarly, an initial 
and a final velocity profile R(t\, r) and R(t2,r), or one den- 
sity profile and one velocity profie determine an LT model. 
Amongst the other options presented, one is to specify a sin- 
gle density (or velocity) profile at a certain time and require 
a simultaneous bang time, which has the effect of remov- 
ing decaying modes. If one is seeking models that are free 
of shell crossings, this must be checked once the arbitrary 
functions have been calculated. 

In this work, the procedure was adapted to the case 
of non-zero A, which means that the analytic evolution ex- 
pressions of the zero-A case had to be replaced by numerical 
determinations. 

We started by choosing both a density profile and a ve- 
locity profile at the present day, solving for the arbitrary 
functions, and then looking at how the model evolved, es- 
pecially the period from last scattering to the present day 
(section 14. H . We found that this procedure typically gave 
severe shell crossings during the evolution, often not far in 
the past, which made the model useless. 

We then tried specifying only the present day velocity 
profile and requiring a simultaneous bang time is = con- 
stant; and we also tried combining a present day density 
profile with a simultaneous bang time. We found this last 
method was particularly good at providing models free of 
shell crossings fsection l4.2H . 

Having found a well-behaved model, it was then possible 
to tinker with it in order to adjust the present day profiles, 
whilst preserving a reasonable evolution (so the bang time 
was often not simultaneous in the final model). The impor- 
tant point is that by considering how physically reasonable 
the evolution of the model is, we can put definite limits on 



the possible present day velocity profiles for a given present 
day density profile, and vice-versa. 



4 RESULTS 

4.1 Present day velocity and density profiles 

Astronomical observations based on different approaches 
provide us with different velocity and density profiles. In this 
section we aim to test the consistency of the various profiles 
with each other and with evolutionary scenarios. Although 
the LT model in principle allows us to match any initial 
velocity profile with any final density distribution, not all 
resulting evolutions are equally natural, and not all are free 
of shell crossings. The velocity and density profiles used be- 
low are those of Figs. |S| and [7] 

• Model 1 — HVP + HM: 

The evolution of this model is presented in Fig. |H| which 
gives the relative density at each point, p(t,r)/pt(t). The 
curves represent the evolution of the derived model without 
a cosmological constant, and the points represent the derived 
model evolution with £7 a = 0.73. 

Because the age of the Universe is different in various 
cosmological models, to compare models with and without 
A, we trace the evolution not in terms of time but in terms 
of redshift. For example, at redshift z = 1 the age of the 
Universe is 5.9 Gy in a A-dust (A-CDM) model, and 4.6 Gy 
in a dust model without A. 

• Model 2 — LVP + HM: 

The evolution results are given in Fig.|S] This model behaves 
differently from what might be expected. The density at 
the center of the SC increases as we trace the evolution 
backwards. This is due to a large central value of the bang 
time function, is, which at the origin is almost 3.5 billion 
years more recent than elsewhere. 

• Model 3 — HVP + LM: 

The evolution of this model is similar to that of model I. 

• Model 4 — LVP + LM: 

This model's evolution is quite like that of model 2. 

The above results suggest that there is no significant dif- 
ference between models with and without cosmological con- 
stant; in other words the value of A is not important on this 
length scale even over cosmological timescales. Additionally, 
models 2 and 4 have the SC growing towards the past, which 
is in contradiction with CMB observations. Therefore there 
is a strong contradiction between the assumption of a low 
amplitude velocity proflie (LVP) and the requirement of a 
reasonable evolutionary scenario. 



4.2 Simultaneous bang time 

The results of the previous section imply that changes in the 
present day velocity profile have a large impact on the evo- 
lution of a model. Moreover, the measurements of peculiar 
velocity have large errors, both observational and systematic 
(we are able to measure only the radial component of the 
velocity). Thus, it is more sensible to obtain a calculated 
velocity profile from other considerations, and compare it 
with observations. 

Measurements of the cosmic microwave background 
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Figure 8. The density evolution of model 1. The centre of the SC 
is at the left, the radius of the LG is at the right, and the radius 
of the GA is near relative distance 0.7. Lines represent models 
with A = 0, points represent models with Q\ = 0.73. 
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Figure 9. The density evolution of model 2. Lines represent mod- 
els with A = 0, points represent models with f^ = 0.73. 



The derived profile of the present day peculiar velocity is 
presented in Fig 1111 As can be seen, the predicted ampli- 
tude of the peculiar velocity around the Shapley concentra- 
tion is around 800 km s" 1 , which is nicely within the range 
between the velocity profiles LVP and HVP. However, this 
model predicts a significant backside infall onto the great at- 
tractor (~ —590 km/s in model the with A = and ~ —620 
km/s in the model with A/0), and very small values for the 
maximum nearside infall onto the GA (ss 235 km/s in the 
A = model, and ~ 250 km/s in the A/0 model). Again 
the curves presented in Fig. 1111 are for models with no cos- 
mological constant, and differences in the peculiar velocity 
profiles between models with and without A are nowhere 
larger than ~ 10%. 

• Model 6 — t B = + LM: 

Fig. llOl shows the evolution of this model, and Fig. llll shows 
the derived present day peculiar velocity profile. Within a 
small distance from Shapley Concentration, the peculiar ve- 
locity is comparable with the LVP, but at all other distances 
the peculiar velocity drops well below observed values. This 
result suggests that the low mass profile (LM) does not pro- 
vide a realistic description of the density distribution around 
the Shapley Concentration. 

• Model 7 — t B = + ~HVP: 

This model aims to recover the high amplitude of peculiar 
velocities around the SC by adjusting the density profile. 
The resulting density profile is obtained from HM by in- 
creasing the density in the SC region. To obtain the peculiar 
velocity maximum near the SC of ~ 1600 km/s, the mass 
of the SC must be almost doubled, i.e. within a radius of 
14.03 Mpc the mass should be 8.5 x 10 15 M , and within 
19.6 Mpc, 19.5 x 10 15 Mq. Such a large mass contradicts the 
astronomical observations mentioned above. 

Therefore, the results obtained in this and the previ- 
ous section suggest a high mass profile (HM) around the 
Shapely Concentration, with peculiar velocity values of 
around 800 km/s. However, all models so far considered, im- 
ply a significant backside infall onto the GA: ~ -700 km/s 
(model 6), ~ -600 km/s (model 5), or ~ -400 km/s (model 
7), which is not satisfactory. 



(CMB) imply that the Universe was very close to homo- 
geneous at the time of last scattering, and as a consequence 
the variation in the bang time function cannot be larger 
than a few thousand years, which in comparison with the 
present age of the Universe is negligible. Therefore, in this 
section the LT model is specified from a present day density 
distribution and the assumption t B — 0. 



4-2.1 Peculiar velocities around the SC 

• Model 5 — t B = + HM: 
The evolution of model 5 is presented in Fig. 1101 For clarity, 
Fig. 1101 only presents the evolution without a cosmological 
constant; evolution with a cosmological constant, as in the 
previous section, does not differ significantly from that with- 
out. Since models with a simultaneous bang time are known 
to describe the evolution of growing modes only, the evolu- 
tion of such models is pretty much as expected, and is free 
of shell crossings. 



4.3 Backside infall 

In this section the issue of backside infall is investigated. As 
stated, there is no observational evidence for any backside 
infall onto the GA. Therefore, models 8 and 9, are intended 
to obtain a backside infall within the range —100 km/s to 
+ 100 km/s. In both models we adopt a density distribution 
around the SC quite similar to the HM profile. 

• Model 8 — t B = + backside infall = -100 km/s: 

To obtain a backside infall onto the GA of —100 km/s, the 
density between the SC and the GA must be around p — 
0.79pi, in the model with A = and around p = 0.78p;, in 
the model with A^0, but otherwise similar to HM. 

• Model 9 — t B - + backside infall = +100 km/s: 

To obtain a backside infall of +100 km/s, the density be- 
tween the SC and the GA must be about p — 0.96pb with 
A = 0, and about p — 0.95pt with A =^ 0, but otherwise 
similar to HM. 

• Model 10 - HM + backside infall = km/s: 

This model aims to check the validity of the results of 
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Figure 10. The evolution of models 5 (labelled "HM") and 6 
(labelled "LM"), with A = 0. The centre of the SC is at the left, 
the radius of the LG is at the right, and the radius of the GA is 
near relative distance 0.7. 
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Figure 11. The present day peculiar velocity profiles used in 
models 5, 6 and 7. The centre of the SC is at the right, the LG is 
at the left, and the GA is near 4000 km/s. 



models 8 and 9, by investigating whether realistic models 
without backside infall really require the density between 
the SC and the GA to be quite close to background value, 
or lower intervening densities are also possible. 

This model was specified by the HM density distribution 
and by a velocity distribution without backside infall (simi- 
lar to that for model 11, see Fig. 1131 . 

If the bang time function ts (r) of this model had turned 
out to have only small variations (e.g. 1000 years), we would 
have found a realistic model without backside infall and with 
low density between the SC and the GA; one whose evolution 
after last scattering would be well behaved. The bang time 
function we actually obtained from this model is presented 
in Fig. 1121 The variation of ts between the SC and the GA 
is almost 700 million years. Since this is far larger than the 
time of last scattering, it results in extensive shell crossings 
and large density variations at last scattering. 

The above results show that the absence of backside infall 




60 80 100 120 140 
Distance R [Mpc today] 

Figure 12. The bang time function tg obtained with model 10, 
against the present day separation of the corresponding world- 
lines. The SC is on the left, and the LG is on the right. 



does indeed require the density between the SC and the GA 
to be close to the background density. 

Models 8, 9 and 10 show that the extent of backside 
infall onto the GA is very sensitive to the density between 
the SC and the GA. If there is no backside infall, the density 
in this region must be p fa 0.85pf,. Models without backside 
infall and with a void between the SC and the GA have 
strong variations in the bang time function, which are not 
consistent with CMB observations. 



4.4 The motion of the Local Group 

So far all our models assumed that the mass inside the sphere 
of radius tlg is the same as the mass in a uniform sphere of 
the same size with a density equal to the background value. 
This assumption, combined with the ts — assumption, 
guarantees that there is no peculiar motion of the Local 
Group towards the GA. Now varying M/E 3 ' 2 at tlg deter- 
mines how much our worldline is gravitationally bound to 
the interior mass, while varying ts only affects how much of 
the evolution determined by M/E 3 ' 2 has happened so far. 
Thus we keep is constant and manipulate M. This induces a 
change of the mass within the sphere of radius tlg, in com- 
parison with the homogeneous background sphere. Thus the 
matching to background happens somewhere further out. 

The assumed density distribution is the similar to HM, 
but instead of a void between the SC and the GA, it has 
density 0.9p;,, so it does not produce any backside infall. We 
then modified this density distribution, as explained below, 
so as to reproduce the LG's 200 km/s velocity toward the 
GA. 

From eq. © it follows that to fit the value of 200 km/s 
for the LG's motion we have to increase M(tlg)- This in- 
crease must be about AM « 95 x 10 15 M if tt A = 0.73, 
and AM « 115 x 10 15 M if fi A = 0. The mass of the homo- 
geneous background within that radius is ~ 1115 x 1O 15 M0 
so the increase is around « 10%. 

This extra mass cannot be added to the SC, i.e. inside a 
sphere of radius 30 Mpc. This is because the SC mass is far 
smaller, and such an increase would cause a larger than ob- 
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Figure 13. The peculiar velocity profiles that result from fitting 
the LG's motion towards GA, using models 11 and 12. The curves 
are labelled by their model numbers. The LG is on the left, and 
the SC is on the right. 



served amplitude in the peculiar velocity profile. Similarly, 
adding the mass in the low density between the SC and the 
GA would lead to an increased peculiar velocity in that re- 
gion of more than 100 km/s, as noted previously. Therefore, 
the only possible way to increase the mass without incon- 
sistency, is to add it at distances comparable with the GA's 
radius or larger. Let us consider these possibilities: 

• Model 11 — increased GA mass: 

A third possibility is to increase the density of shells at the 
radius of the GA. The density contrast at the GA's radius 
was increased from the initial 6 — 1 to 6 — 1.4, thus induc- 
ing a change of the Great Attractor mass 2 from an initial 
M « 4.7 x 10 15 M Q to M « 5.9 x 1O 15 M . The velocity dis- 
tribution around the GA predicted by this model is in quite 
good agreement with the HVP (but only around the GA). 

• Model 12 — increased LG mass: 

The last possibility is to increase the density of the region at 
the radius of the Local Group by AS ~ 0.1. The prediction 
of this model for the velocity distribution around the GA is 
also in quite good agreement with the LVP (but only around 
the GA). 

The deduced peculiar velocities for these models are 
presented in Fig. 1131 where for clarity, only results from the 
A/0 model are plotted. 

As can been seen the peculiar velocity around the GA 
is very sensitive to its mass. A mass change of order 25% 
leads to a peculiar velocity change of about 50%. So if we 
could extract from observation the exact value of the pecu- 
liar velocity profile around GA, we could estimate its mass 
to good precision. However, due to the GA being in the Zone 
of Avoidance, this is a very difficult task. 



2 To estimate this mass we take a roughly cubic blob out of the 
spherical density maximum near radius 130 Mpc from the SC. 



4-4-1 The SC's contribution to the LG motion 

The motion of the Local Group with respect to the CMB is 
around 600 km s _1 . The cause of this motion is still a riddle. 
The contribution to this motion from different superclusters 
and the Local Void have been systematical studied. How- 
ever, there still remains the question of the Shapley concen- 
tration's contribution to this motion. For example, Tonry et 
al. (2000) suggest that the SC's effect is rather negligible. 
On the other hand X-ray observations (Kocevski, Mullis & 
Ebeling 2004) imply that the SC makes a significant contri- 
bution to the motion of the Local Group. 

To investigate this issue, we considered a model in which 
the density distribution is that of the Shapley Concentration 
in the HM profile out to where its relative density drops to 
1 (around 30 Mpc), and then remains at the homogeneous 
background value beyond that. The peculiar velocity thus 
obtained is presented in Fig. 1141 It is apparent that the mo- 
tion towards the SC at the distace of the LG is very low, 
approximately 35 km/s. 

The explanation for this phenomenon can be illustrated 
by the following approximate calculation based on the LT 
evolution formula J3J. Let us consider a uniform gravitating 
sphere of radius R and mass M, with a small extra mass m 
(e.g. the SC) at the center, m « M. The influence of this 
additional mass on the velocity of test particles at distance 
R is estimated as follows: 



-AV = - (V - V b ) = 
c c 



2£+ 2(M + ™) + l A7?2 + 



2M 1 . „„ 
-W + 3 AR 



2E- 



,„„ 2M l (mu , 2mc 2 . 

( 2S + 1T + 3 A * 2)(1 + ^F ) + 



2E+2 -f + l AR2 



RR 



R 2 H 



(6) 



In the above formula, a negative velocity means motion to- 
wards the center, hence the peculiar velocity at radius R is 
towards the extra mass. Converting the above relation to a 
ready-to-use formula, we obtain: 

M 



AV = 4.3 x 10 4 



D 2 h' 



(7) 



where M is the extra mass in units of 10 Mg, D is the 
distance in Mpc, h is defined by Ho — 100 h km/s/Mpc, and 
AV is the amount of infall velocity due to the central mass 
excess in km/s. Therefore, although the mass of the Shapley 
Concentration is large (M ~ 10), the distance to it is also 
large, (D ~ 190) and hence its contribution to the motion 
of the Local Group is small. 



5 CONCLUSIONS 

We have developed a relativistic approach to modeling the 
Great Attractor and the Shapley Concentration, with the 
best models being models 11 and 12. The Great Attractor 
region has been the subject of many observational and the- 
oretical studies. However, due to the Great Attractor being 
in the Zone of Avoidance, all of these studies had to make 
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Figure 14. The peculiar velocity profile for a Universe consisting 
of the SC added to an otherwise uniform density. The SC is at 
the right. 



some assumptions. The POTENT analysis (Bertschinger & 
Dekel 1989; Dekel, Bertschinger & Faber 1990a, 1990b; Ko- 
latt, Dekel & Lahav 1995) assumed Newtonian gravity and 
no vorticity. Tonry's infall analysis (Tonry et al. 2000) as- 
sumed a spherical infall onto the Great Attractor, which 
implied the existence of backside infall. 

Our approach instead assumes spherical symmetry 
around the Shapley Concentration. This assumption allows 
us to investigate whether or not a flow towards Great At- 
tractor can be a part of a larger flow towards the Shapley 
Concentraton. Although, our spherical assumption is only a 
first approximation, it is not seriously off, especially for the 
region this side of the SC. Our research should be viewed 
as one of several complementary approaches to studying the 
Great Attractor. The overlap of all these approaches should 
provide us with better understanding the Great Attractor 
region. 

The conclusions of our study of the SC, GA and LG 
region are as follows: 

(i) The peculiar velocity around the SC is ~ 800 km/s 
towards its centre. 

(ii) The lack of evidence for backside infall onto the GA 
implies that the density between the GA and the SC must 
be about S « 0.9. 

(iii) The mass of the GA is around 4 - 6 x 10 15 M . 

(iv) The SC's contribution to the LG peculiar motion is 
negligible. 

(v) On the scales of this investigation, the value of the 
cosmological constant does not have a significant impact on 
the formation and evolution of structure. 
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